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ABSTRACT 

We present spectroscopy and photometry of GD 448, a detached white dwarf- M dwarf 
binary with a period of 2'M7. We find that the NaI8200A feature is composed of 
narrow emission Unes due to irradiation of the M dwarf by the white dwarf within 
broad absorption Hnes that are essentially unaffected by heating. Combined with an 
improved spectroscopic orbit and gravitational red shift measurement from spectra of 
the Ha line, we are able to derive masses for the white dwarf and M dwarf directly 
(0.41 ± O.OlMoand 0.096 ± 0.0041/©, respectively). We use a simple model of the 
Ca II emission lines to establish the radius of the M dwarf assuming the emission from 
its surface to be proportional to the incident flux per unit area from the white dwarf. 
The radius derived is 0.125 ± 0.020 i?© . The M dwarf appears to be a normal main- 
sequence star in terms of its mass and radius and is less than half the size of its Roche 
lobe. The thermal timescale of the M dwarf is much longer than the cooling age of the 
white dwarf so we conclude that the M dwarf was unaffected by the common-envelope 
phase. The anomalous width of the Ha emission from the M dwarf remains to be 
explained, but the strengh of the line may be due to X-ray heating of the M dwarf 
due to accretion onto the white dwarf from the M dwarf wind. 

Key words: white dwarfs - binaries: spectroscopic ~ binaries: close - novae, cata- 
clysmic variables - stars: individual: GD448 - stars:low-mass, brown dwarfs 



1 INTRODUCTION 

GD448 (= WD 0710+741, LP034-185) was shown by 
Marsh & Duck (1996; MD96 hereafter) to be a detached 
white dwarf- M dwarf binary with a period of 2M7. This 
is the shortest period known for such a binary and places 
it in the centre of the "period gap" from 2 to 3 hours in 
which very few cataclysmic variable stars (semi-detached 
white dwarf- main-sequence binaries) are found. From the 
cooling age of the white dwarf, MD96 found that GD 448 
was born in the period gap and had never been a cata- 
clysmic variable star. They also found that the Ha emission 
line from the heated face of the M dwarf is much broader 
than expected. Since the thermal timescale of the M dwarf is 
much longer than the cooling-age of the white dwarf the ra- 
dius of the M dwarf in GD 448 is essentially unchanged from 
the radius at the end of the common-envelope phase. This 
enables us to study the effect of the common-envelope phase 
on the M dwarf. Unfortunately, the data of MD96 were of 
insufficient quality to establish the radius of the M dwarf. 

In this paper we present V and I band lightcurves of 
GD 448 and much improved spectroscopy which has enabled 



us to determine the radius of the M dwarf in GD 448 and to 
improve the mass estimates for both components. 



2 OBSERVATIONS 
2.1 Spectroscopy 

We observed GD448 over four nights (11-14 January 1996) 
using the double-beam spectrograph ISIS on the 4.2m 
William Herschel Telescope. Data were obtained over three 
nights in variable conditions. We used the R1200R grating 
in the blue arm and the R600R grating in the red-arm with 
TEK CCD detectors to obtain spectra covering 400A around 
Ha at a dispersion of 0.4A per pixel and an 800A region cov- 
ering the Na 1 8200 A doublet and Call infrared triplet at a 
dispersion of 0.8A per pixel. A total of 161 spectra were 
obtained in each spectral region. A 0.9" wide slit which 
projected onto 1.6 pixels at the detectors was used to en- 
sure accuracy in the measured radial velocities. Exposure 
times were 300s -500s. Optimal extraction was used to ex- 
tract spectra from our CCD images. All the observation were 
interspersed with observations of a copper-argon arc at least 
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Table 1. Position of GD448 and comparison stars Coordinates 
are equinox 2000.0, epoch 1996.24 and are accurate to a few arc- 
seconds. 

Star RA Dec 

GD448 07''17™09.6'' +74°00'42" 
CI 07''17™03.9= +74°03'32" 

C2 07''16™-41.3= +74°00'23" 

Figure 1. The V and I lightcurves of CD 448 (upper and lower 
panels, resp.). 
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every 40 minutes and before and after every re-pointing of 
the telescope. Spectra of the arc were extracted at the same 
detector position as the object spectra. The wavelength scale 
determined from a polynomial fit to measured arc line posi- 
tions was interpolated in time for the object spectra. 

The spectra were flux calibrated using spectra of 
BD-l-26 2606 (Oke 1983) to correct for the wavelength- 
dependent instrument sensitivity. No correction for slit 
losses was attempted so the spectra have an arbitrary ab- 
solute fiux scale. Telluric features in the red-arm data were 
removed using the technique of Wade & Horne (1988). The 
rapidly rotating B-type star HD 84937 was used to determine 
the telluric absorption spectrum. The removal is naturally 
imperfect as the strength of the telluric absorption had to 
be estimated from the airmass at which each spectrum of 
GD 448 was observed. 



2.2 Photometry 

We used the Im Jacobus Kapteyn telescope to obtain V and 
I band CCD images of GD 448. The detector used was a TEK 
CCD with 1024^ pixels giving an image scale of 0.34" per 
pixel. Observations were obtained during the first 1-2 hours 
of the nights 28 March to 2 April 1996. A total of 84 useful 
images were obtained with the I filter and 30 in the V filter. 
Normalized average images of the twilight sky were used to 
apply fiat-field corrections to all the images. 

We used profile fitting to determine differential magni- 



Figure 2. The measured radial velocities for the Call triplet 
emission lines (upper panel) and adopted fit. 
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tudes of 8 stars in the images including GD 448, all of which 
were used to determine the shape of the profile to be used in 
each image. After careful inspection of the images and the 
resulting magnitudes of these stars we chose the comparison 
star CI and a check star C2 (= CSC 4372 00331) whose po- 
sitions are given in table |l[ The positions of these stars and 
GD 448 were measured from our images using positions for 
5 Hubble Guide Star Catalogue stars to establish the lin- 
ear astrometric solution. The magnitude difference between 
these stars is seen to be constant to within « 0™01. The 
lightcurve of GD 448 with respect to CI is shown in Fig. ^ 
where the phase has been calculated using the ephemeris 
given below. The I band lightcurve clearly shows the effect 
of the irradiation of the M dwarf by the white dwarf. The 
incomplete coverage of the V lightcurve makes it difficult to 
judge whether the variabilty at this wavelength is due to ir- 
radiation, although a cosine function gives a good fit to the 
incomplete lightcurve and suggests the amplitude of the V 
lightcurve is half that of the I lightcurve. 



3 RESULTS 

3.1 The Call triplet 

The red-arm spectrum shows the Call triplet emission lines 
seen by MD96. These are of roughly equal strength and 
show sinusoidal variations of intensity and Doppler shift 
with phase. This can be seen in Figs. |2|, ^ and ^ which 
are described in more detail below. 

We used a series of Gaussian fits to individual spectra 
to measure the radial velocity of the Ca II triplet in a similar 
manner to MD96. The spectra are normalised to a contin- 
uum value near the Call lines determined from a cosine 
fit to the I-band lightcurve. We find Ci = 1 — 0.03 cos((/>), 
where Ci is the continuum value and (j) is the orbital phase 
in radians, phase zero being times at which the M-star is 
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Figure 3. Trailed and phase-binned spectra of the Call 8542A Une (left panel), the model spectra for the simple model of Call emission 
and the adopted value / (middle panel) and the H« emission line after subtraction of the white dwarf absorption (right panel) . 
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Table 2. The parameters of the Gaussian profiles used to fit the 
Call emission lines. 



parameter 
FWHM(kms-i) 

h(8662)/h(8542) 
h(8498)/h(8542) 
Ai 



Value 
58.6 ±0.5 

0.830 ±0.010 
0.708 ±0.007 
0.409 ± 0.002 
-0.267 ±0.003 



Notes 

Mean width of lines. 
No phase variation seen. 
Ratio of heights 
Ratio of heights 



1: Fit to height dependence on phase (f>, h = A + i?cos(0) 



Table 3. Circular orbit fits 

Spectral line(s) 7(kms~"'^) ii'(kms~"'^) 

Call IR triplet emission -40.8 ± 0.35 123.1 ± 0.4 

Ha absorption -24.2 ± 0.6 -33.5 ±0.9 

Ha emission -38.1 ± 0.4 125.4 ±0.6 

Na 1 8200A absorption -38.6 ± 1.5 142 ± 3 

Na 1 8200A emission " 120 ± 3 

88OO.3A absorption — 148 ± 5 



closest to the Earth. The ratio of the heights of the Call 
lines was held fixed at a value determined from a previous 
fit and the variation of the height of the lines with phase was 
set by a fit to the measured individual heights of the form 
A + Bcos{(j)). We detected no variation of the width of the 
lines with phase so this height variation has the same form 
as the equivalent width variation seen in Fig. i. Values of the 
parameters used in the fit are given in table 4 The reduced 
chi-squared values for the fits were typically 0.9-1.1. 

We used a periodgram of our radial velocity measure- 
ments combined with those of MD96 to confirm the orbital 



period, P. The cycle count between the two data sets is 
confirmed by the lightcurve. From the sine curve fit of the 
form 7m + -R'm^^^ sm{(j)) to the combined radial velocity data 
shown in Fig. ^ we find the following ephemeris: 

HJD(ro) = 2450094.69006(7) + 0.10306437(7)£, 

where figures in parentheses are l-cr uncertainties in the final 
digit and To refers to dates at which the M-star is closest 
to the Earth. The values of 7m and i^M^^' determined from 
this fit are given in table |[ This ephemeris agrees with that 
of MD96 once a superfluous digit in their quoted period is 
accounted for. 

The average red-arm spectrum after subtraction of the 
Doppler shift due to this orbit is shown in Fig. ^. In addition 
to the Call triplet, there are several other emission lines vis- 
ible. Some of these are labelled with tentative identifications 
which are based on a simple model of emission from an op- 
tically thin gas with a range of effective temperature around 
5000- lOOOOK. Also visible is the NaI8200A feature, which 
is distorted by the averaging process, as well TiO and other 
absorption features. 



3.2 The Na I doublet. 

We experimented with several models to account for the 
variation with phase of the NaT doublet. The most successful 
and easiest to interpret is that of narrow emission lines, sim- 
ilar to the Call emission lines, with semi-amplitude K^'^^ 
superimposed on broad absorption lines which are constant 
with phase, i.e., with the same semi-amplitude as the centre- 
of-mass of the M dwarf, Km- To measure the values of K^^^ 
and Km we used a fitting method similar to Moran et al. 
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Figure 4. The average red-arm spectrum after removal of the Call radial velocity Doppler shift. 
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Table 4. The parameters of the Gaussian profiles used to fit the Nal doublet. Parameters in bold type are fixed quantities. 

Wavelength(A) 8183.3 8194.8 8183.3 8194.8 8207.5 ± 0.4 

FWHM(A)i 3.3 ±0.2 3.8 ±0.2 1.8 ±0.1 2.5 ±0.9 

Height -0.079 ±0.004 -0.073 ± 0.004 1.07 ±0.09 0.99 ± 0.08 -0.010 ± 0.002 

A — — 0.066 ±0.003 0.057 ±0.003 — 

B — — -0.038 ±0.004 -0.048 ±0.004 — 



(1997), in which a simultaneous Gaussian fit to all the spec- 
tra is used to determine the orbital parameters. A single 
Gaussian was used to model each emission and absorption 
component for both Nal lines and a third Gaussian was used 
to account for weak absorption line nearby. The parameters 
of these Gaussians are given in table ^. The variation in 
height of the emission lines was established from a cosine fit 
to measured individual heights, note that this is similar to 
the variation seen in the Call emission lines. The reduced 
chi-squared value of the final fit is 1.05 (9285 data points) 
which is acceptable given the difficulties in removing telluric 
features in this region. The observed NaT line spectra and 
fits have been coadded into phase bins to produce the trailed 
spectra shown in Fig. ^ . 

We were able to confirm the value of T^m using the ab- 
sorption line seen near 8800A in Fig. M, which we tentatively 
identify as Fel 8801.8 A. We used a fit of a single Gaussian 
to this line to determine Km- The line shows no variation 
of depth or width with phase and is unaffected by telluric 
features. The nearby Mgl emission line was accounted for in 
the fit and was found to behave in a similar way to the Ca II 
lines. We find Km. = 148±5kms~^, in good agreement with 
the value determined from Na I. The weighted average of the 
two values is 144 ± 3km s~^, which is our adopted value for 
Km. 



Table 5. The parameters of the Gaussian profiles used to fit H«. 

FWHM(A) 
2.08 ±0.02 
1.321 ±0.034 



Component 
Emission line 
Absorption component 1 
Absorption component 2 
Absorption component 3 
Absorption component 4 
Absorption component 5 
a: The actual height of the 
multiplied by the cosine 



3.3 The Hq line. 



6.87 ±0.30 
28.36 ±0.90 
67.35 ±2.41 
153.6 ±4.40 



Height 

1.030 ± 0.0006" 
-0.175 ±0.002 
-0.070 ± 0.002 
-0.109 ±0.003 
-0.117 ±0.003 
-0.074 ± 0.003 



emission is given by this value 
variation described in the text. 



The contribution from the M dwarf to the continuum near 
Ha is negligible so the continuum was set to unity using 
a linear fit to the continuum 8000-8500km s~^ either side of 
the line. We used five Gaussian profiles to model the ab- 
sorption profile and a single Gaussian profile to measure the 
emission feature. A series of fits were used to establish the 
height variation of the emission line, which is similar to the 
Call emission. We also established the widths and heights 
of the Gaussian profiles given in table and the apparent 
systemic radial velocity, 7m and semi-amplitude, K^'^ of 
the Ha emission line given in table ^ using a simultaneous 
fit to all the spectra. 
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Figure 5. Trailed and pliase-binned spectra around tiie Nal 
8200A absorption line (left panel) and the best fit using multiple 
Gaussian fits (right panel). 
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where Mwd is the mass of the white dwarf, Mm is the mass 
of the M dwarf and i is the incHnation. The actual mass of 
the white dwarf can be estabhshed from the gravitational 
red shift and the mass -radius relation for white dwarfs. To 
the difference (7wd — 7m) = 13.9 ± 0.6km s~^ measured 
from the Ha line must be added corrections for the red 
shift of the M dwarf (0.5km s~^) and the difference in trans- 
verse Doppler shifts (0.1km s~^) as described by MD96, as 
well as for the potential at the M dwarf due to the white 
dwarf (0.35km s~^) and vice versa (-0.08kms~i), a correc- 
tion which MD96 omitted. The mass of the white dwarf 
is then determined from the model cooling curves for low- 
mass helium white dwarfs of Althaus & Benvenuto (1997). 
For white dwarfs with an effective temperature of 19 000K 
and a mass near 0.4Mq , we find the mass in solar masses 
is related to the gravitational red shift, Vg, in kms~^ by 
AI — 0.181 -I- Q.0157f5. The inclination can then be deter- 
mined using equation (1) of MD96. The corrections to the 
gravitational red shift depend slightly on the assumed val- 
ues of i and so an iterative scheme is used to find consistent 
values for all the parameters. The dependence of the results 
on the assumed radius of the M dwarf is negligible for any 
reasonable estimate. The inclination and masses derived are 
given in table The radius of the white dwarf is also given, 
although no error is quoted as this is strongly correlated 
with the uncertainty in the mass. 



Table 6. Summary of the parameters of GD448. 

Parameter white dwarf M dwarf 

Mass {Mq ) 0.41 ± 0.01 0.096 ± 0.004 

Radius(_RQ) 0.018 0.125 ±0.02^0 

Inclination(°) 29.3 ± 0.7 

Period (d) 0.10306437 



We were concerned that the variation of the Ha emis- 
sion line might affect our estimates of 7wd and Kwd - the 
apparent systemic radial velocity and semi-amplitude of the 
white dwarf. To avoid this we excluded a region on either 
side of the predicted emission line position from the fitting 
process with similar width to the FWHM of the emission 
line (~ 100km s~^). We then used a simultaneous fit to all 
the spectra to measure the 7wd and Kwd- The results vary 
slightly depending on the width of the excluded region, and 
this additional uncertainty has been included in the adopted 
values shown in table ^ 

Our fit to the white dwarf absorption line was used 
to subtract the contribution of the white dwarf from the 
Ha line. The remaining M dwarf emission line spectra were 
coadded in phase-bins to produce the trailed spectra shown 
in Fig. bl The equivalent width variation of this line is shown 
in Fig.^ It is obvious from Fig. ^ that the Ha emission is 
much wider than the CaH emission lines, a feature noted 
by MD96. Less obvious but still visible is the double-peaked 
shape of the Ha emission near phases 0.25 and 0.75. These 
anomalies are discussed in more detail in section ^ 

3.4 The mass of the white dwarf and the M dwarf. 

The measured values of Kwo, Km and the orbital period 
lead directly via Kepler's laws to AfwD sin'' i and Mm sin^ i. 



4 DETERMINATION OF THE RADIUS OF 

THE M DWARF FROM THE CA II TRIPLET. 

The CaH emission in GD448 is due to emission from the 
heated face of the M dwarf so K^^^^ is lower than Km- The 
correction AK = Km — ^m''^' depends on the pattern of the 
emission over the surface of the M dwarf and its radiative 
properties (e.g., optical depth, limb darkening) as well as 
the geometry of the system. The only free parameter in the 
system geometry is the filling factor of the M dwarf, /, the 
ratio of the radii of the M dwarf and its Roche lobe measured 
from the centre-of-mass to the inner Lagrangian point. A 
value of / = 1 corresponds to a Roche lobe filling star. The 
other parameters (e.g .. inclination, mass ratio) have been 



derived in section 3.4 above. 

For our model of the emission from the M dwarf, we 
assumed a source function of the form 1 + er^ where r is 
the optical depth along the line of sight and jj, is the cosine 
of the angle between our line of sight and the normal to 
the surface. The dependence of the emission line flux from a 
point on the star, F, with viewing angle then has the form: 

F oc (1 + e ■ m)(1 - 6""°/") - e ■ roe"""/'', 

where to is the vertical optical depth through the emitting 
region, which is assumed to be constant over the surface of 
the star. In the optically thick case (to ^1), e is related to 
the standard linear limb darkening parameter which is given 
by e/(l + e) when e is positive. In cases where e is negative, 
i.e., source function increasing with height, as might be ex- 
pected in the case of an irradiated atmosphere, tq must be 
less than — 1/e to avoid negative source function values. 

We used numerical integration over a model star de- 
fined by a surface of constant Roche potential to predict 
emission line strengths and radial velocities for various val- 
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Figure 6. Variation of the equivalent width of the Ho (filled circles) and CaII8542A (crosses) emission lines with phase. Solid lines 
are fits of the form A + B sin((/)) + C sin(2(/)). The dotted line is the best fit for an optically thick model to the Call 8542A line. 
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ues of /, e and ro assuming the emission line strength to 
be proportional to the incident flux per unit area from the 
white dwarf. These were calculated at the same phases as 
the observed data and given the same weighting in a sine fit 
that was used to calculate /S.K. This has an almost linear 
dependence on / and so it is very easy to find / which gives 
the correct value of /S.K. 

A further constraint on our model of the Ca II emission 
comes from the lightcurve of the Call lines, i.e., the variation 
of equivalent width with phase. By fitting a function of the 
form A + i3sin(0) + Csin(2(/)) to the observed equivalent 
width variation, we find C/B = -0.02 ± 0.02. This ratio 
describes the magnitude of any non-sinusoidal component 
of the lightcurve (i.e., its shape) independent of any zero- 
point or scaling applied. For each value of e and ro used, the 
predicted equivalent widths for the optimum filling factor 
/o were calculated at the same phases as the observed data 
and given the same weighting in a fit of the same function. 
The ratio C/B from these fits rules out models with values 
of To ^ 1.0. These models show a non-sinusoidal component 
that is much larger than that observed {C/B ~ 0.1), i.e., 
they are the wrong shape. This is shown in Fig. ^ 

The values of /o found for valid values of ro vary from 
0.46 (ro < 1) to 0.49 (ro = 1) whatever the value of e. 
The uncertainty in Km leads to an further uncertainty of 
0.06 in /o. The absolute radius of the M dwarf is then 
0.125 ± 0.015 ± 0.005 7?0 , where the uncertainties quoted 
are the random error inherited from Km and the systematic 
error due to the range of /o values, respectively. The com- 
bination of systematic and random uncertainties is a moot 
point, but the result cannot be worse than a simple linear 
addition and so we will adopt a value of 0.125 ± 0.02 i?© . 



The relative sizes of the white dwarf, the M dwarf and the 
Roche lobe of the M dwarf are shown in Fig. ^ 

The position of the M dwarf in the mass-radius plane 
is shown in Fig. ^ Also shown in Fig. ^ for comparison are 
masses and radii for other M dwarfs given by Clemens et al. 
(1997). Their masses are based on their revised empirical 
color-luminosity relation and the luminosity-mass relation 
of Henry & McCarthy (1993) based on 37 visual binaries. The 
radii are derived from their updated colour-temperature and 
colour-bolometric correction relations based on a volume 
limited sample of 127 M dwarfs. The agreement between 
their indirect method and our more direct method is quite 
satisifactory. 



5 DISCUSSION. 

Comparing the position of CD 448B in the mass-radius 
plane to those of other M dwarfs, we see that it appears to 
be a perfectly normal main-sequence M dwarf. The thermal 
timescale for the M dwarf (3 x lO^y) is much longer than the 
cooling age of the white dwarf (5 x lO^y) , so the radius will 
not have changed significantly since the common-envelope 
phase. The argument of MD96 which suggests that CD 448 
was born in the period gap, i.e., has never been a CV, is not 
affected by anything presented here. 

The anomalous width of the Ha emission line is partic- 
ularly strange given that its intensity variation with phase 
and radial velocity semi-amplitude are almost identical to 
the much narrower Call emission lines. For example, from 
a fit of the form A -I- Bsin{(j>) + C sin(20) to the equiva- 
lent width variation of the Ca II 8542A line shown in Fig. ^ 
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Figure 7. A cross section tlirougli the orbital plane of GD 448 
showing the relative sizes of the white dwarf (filled circle) and the 
M dwarf (partially hatched) with its Roche lobe (dashed lines). 
The hatched region of the M dwarf shows approximately where 
the Call and Ha arise. The + symbol marks the centre of mass of 
the system and the X symbol the centre of mass of the M dwarf. 
The axes are in units of the orbital separation (rx and ry ) and the 
amplitude of the sine and cosine components of projected radial 
velocity (K^ and Ky). 
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Figure 8. GD 448 in the mass-radius plane (filled circle) com- 
pared to other low-mass stars and brown dwarfs from Clemens 
et al. (open circles). 
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we find, BjA^ -0.68 ± 0.01 and C/B = -0.03 ± 0.02. The 
other Ca II lines give similar results. For the Hq emission we 
find the remarkably similar values B jA = —0.71 ± 0.02 and 
C jB = —0.05 ± 0.02, i.e., the Call and Ha emission lines 
have nearly identical amplitudes and shapes. This appears 
to rule out any broadening mechanism requiring a distribu- 
tion of Ha emitting material much different from the Ca II 
emitting region, e.g., photospheric emission. 

Our favoured explanation for the shape of the Ha line 



is some broadening mechanism (e.g., natural damping or 
Stark broadening) within an optically thick medium. This 
is at variance with the results of our simple model for the 
emission of Ca II which appears to rule out emission from an 
optically thick gas. However, our simple model is not com- 
plete as it does not account for any variation of the optical 
depth over the surface of the M dwarf, nor for the finite 
geometrical depth of the emitting layer. In the case of emis- 
sion from an optically thin gas, the emission line strengths 
are expected to be proportional to the oscillator strengths 
of the transitions. In the case of the Call lines, these ratios 
are Ca II 8498 : Ca II 8542 : Ca II 8662 = 1:9:5. The observed 
ratio of the emission line strengths is 1 : 1.4 : 1.2, which sug- 
gest that at least some part of that emission is due to op- 
tically thick material. The double-peaked shape of the Ha 
emission might then be explained as self-absorption of the 
line. Clearly, more detailed modeling would be required to 
demonstrate that these suggestions are feasible, but if they 
do adequately explain the shape of the Ha and Call emis- 
sion lines, these lines would then be giving information on 
the density and temperature structure within the emitting 
region. 

If our assumption that the pattern of Call emission is 
proportional to the irradiating fiux is wrong, it might seem 
possible to explain the width of Ha as rotational broadening 
of an intrinsically narrow line. This would then require a very 
odd distribution of Ha and Call emission over the surface 
M dwarf to explain the shape of the Ha line. In fact, we 
used a very general model of the Ha emission in which the 
pattern of emission over the surface of the M dwarf was 
allowed to vary independently at every grid point, each of 
which emits a narrow line. The model includes the effects of 
orbital phase smearing and the instrumental resolution. We 
then attempted to use a maximum entropy reconstruction 
to determine a pattern of emission that produced emission 
lines consistent with the data, but were unable to find any 
such pattern for any value of /. This is a similar problem 
to that found by Rutten & Dhillon (1994) in the case of the 
cataclysmic variable star DW UMa, although that case is 
complicated by the presence of an accretion disc. 

The photometry for GD448 has only been used to con- 
firm the correct orbital period for the system and to nor- 
malise the continuum value for the red-arm spectra. A sim- 
ple model of the irradiation shows that the amplitude of the 
lightcurve is consistent with that which would be expected 
from an irradiated star, but given the uncertainties in the 
heating mechanism and the small amplitude of the modu- 
lation this gives no constraint on the system geometry in 
practice. 

The problem of the strength of the Ha emission dis- 
cussed by MD96 remains and is made worse by the small 
radius for the M dwarf derived here. Briefly, there is in- 
sufficient radiation below the Lyman limit from the white 
dwarf incident at the M dwarf surface to produce the ob- 
served emission line strength by photoionisation of hydro- 
gen. A similar problem has been identified for the detached 
white dwarf- M dwarf binaries PG1026-I-002 (Saffer et al., 
1993) and WD2256-H249 (Schmidt et al., 1995). MD96 sug- 
gest that chromospheric activity may excite hydrogen atoms 
to the n=2 level, from where they are ionised. 

An alternative explantion comes from the possibility 
of accretion onto the white dwarf from the M dwarf wind. 
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a process seen in the white dwarf- K dwarf detached bi- 
nary V471 Tau (Mullan et al. 1991). The resulting irra- 
diation from the soft X-rays produced would be sufficient 
to produce the observed Ha flux for mass-loss rates only a 
few times greater than the Solar mass-loss rate. To calcu- 
late the Ha flux, we first used the mass - colour - absolute- 
magnitude relations for M dwarfs from Clemens et al. (1998) 
to find (V-I)=4.0 ±0.1 and My = 16.2 ± 0.3. We then used 
Mv =10.15 for the white dwarf (Bergeron et al., 1992) and 
(V-I)=-0.204 for a hydrogen-rich white dwarf of 19 000K 
(Bergeron et al.l995) to predict a contribution of 15%±3% 
from the M dwarf in the I band (~ 8200A). This is in perfect 
agreement with the estimate of MD96. We used (R-I)=2.16 
for the M dwarf (Bessel, 1991), which implies a spectral type 
of M5.5 -M6, and (R-I)=-0.11 for the white dwarf (Bergeron 
et al. 1995) to flnd a contribution of 2.2%±0.7% from the 
M dwarf in the R band (~ 7000A). The peak equivalent 
width of the Hq emission relative to the combined contin- 
uum is I.OA so the true equivalent width of the emission is 
~ 45A. The zero-point of the R-band magnitude scale for 
M dwarfs is approximately 30 x 10"" Wm"^A"^ (AUard & 
Hauschildt. 1995) and the V magnitude of GD 448 is 14.97 
(Bergeron et al., 1992) so the observed Ha flux from the M 
dwarf is ~ 3 X 10~^*Wm~^. For the purposes of this dis- 
cussion we can assume this Ha flux is radiated isotropically 
from the M dwarf, in which case the luminosity in Ha of the 
M dwarf is 3 x 10^° W. If all the X-radiation emitted from the 
white dwarf intercepted by the M dwarf (1.5% for isotropic 
emission) is re-processed as Ha emission, the X-ray luminos- 
ity of the white dwarf is 2 x 10^^ W. The mass accretion rate 
required to provide this energy is ~ 1 x lO~^''M0y~^. The 
efficiency of Bondi-Hoyle accretion from the M dwarf wind is 
expected to be quite high (~ 20%, Mullan et al. 1991). Fur- 
thermore, models of X-ray heating of dwarf M-star chromo- 
spheres suggest that perhaps half of the energy is re-emitted 
as Balmer emission if the X-rays are absorbed high in the 
atmosphere (Cram 1982). The total mass- loss rate required 
to produce Ha could then be as low as ~ 1 x 1O~^^M0 y~^ - 
comparable to the Solar mass-loss rate and much lower than 
the observed mass-loss rate in other M dwarfs (Mullan et al. 
1992). The observed X-ray flux in this scenario would be 
~ 2 X 10"^®Wm"^ = 2 X 10""ergss"^cm"^, which is weU 
within the capabilities of the current generation of X-ray 
satellites. 



REFERENCES 

Allard F., Hauschildt P. H., 1995, ApJ, 445, 433 
Althaus L. G., Benvenuto O. G., 1997, ApJ, 477, 313 
Bergeron P., Saffer R. A., Liebert J. ,1992, ApJ, 394,228 
Bergeron P., Wesemael P., Beauchamp A., 1995, PASP, 107, 1047 
Bessel M. S., 1991, AJ, 101, 662 

Clemens J. C, Reid I. N., Gizis J. E., O'Brien M. S., 1998, ApJ, 

496, 352 
Cram L. E., 1982, ApJ, 253, 768 
Henry T. J., McCarthy D. W. 1993, AJ, 106, 773 
Marsh T. R., Duck S. R, 1996, MNRAS, 278, 565 (MD96). 
Moran C, Marsh T. R., Bragaglia A., 1997, MNRAS, 288, 538 
Mullan D. J., Shipman H. L., MacDonald J., Sion E. M., 1991 

ApJ, 374, 707 

Mullan D. J., Doyle J. G., Redman R. O., Mathioudakis M., 1992, 

ApJ, 397, 225 
Oke J. B., 1983, ApJ, 266, 713. 

Rutten R. G. M.; Dhillon V. S., 1994, A&A, 288, 773 

Saffer R. A.; Wade R. A.; Liebert J., Green R. P., Sion E. M., 

Bechtold J., Foss D., Kidder K. 1993, AJ 105, 1945 
Schmidt G. D., Smith P. S., Harvey D. A., Grauer A. D., 1995, 

AJ 110, 398 
Wade R. A., Home K., 1988, ApJ, 324, 411. 



6 CONCLUSION. 

Extensive analysis of spectra covering the Ha line, the Ca II 
triplet and the Nal doublet from GD448 have enabled us to 
determine the mass and radius of the M dwarf companion 
to the white dwarf GD 448. Its position in the mass-radius 
relation places it squarely on the main-sequence. This sug- 
gests that the common envelope phase had little effect on 
the structure of the M dwarf. The anomalous width of the 
Ha emission from the M dwarf remains to be explained. 
The anomalous strength of the Ha line may be due to X- 
ray heating of the M dwarf due to accretion onto the white 
dwarf from the M dwarf wind. 
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